function __test_natural_spline()
	% functie de test pentru functia de interpolare folosind spline-uri
        % naturale

	%Testul nr.1
	x=[0 1 2 3]
	y=[1 1 3 0]
	a=1.5
	val=natural_spline(a,x,y);
	printf("Valoarea polinomului este:\n");
	disp(val);
	%solutia este:  val= 2.22


	%Testul nr.2
	x=[0 1 2 3]
	y=[-1 0 2 1]
	a=1.5
	val=natural_spline(a,x,y);
	printf("Valoarea polinomului este:\n");
	disp(val);
	%solutia este:  val=1.15


	%Testul nr.3
	x=[1 2 4]
	y=[3 4 6]
	a=1.5
	val=natural_spline(a,x,y);
	printf("Valoarea polinomului este:\n");
	disp(val);
	%solutia este:  val=3.5000


	%Testul nr.4
	x=[1 2 3]
	y=[5 -1 4]
	a=1.5
	val=natural_spline(a,x,y);
	printf("Valoarea polinomului este:\n");
	disp(val);
	%solutia este:  val= 0.96

	%Testul nr.5
	x=[-1 0 1 2]
	y=[-2 -1 2 4]
	a=1.5
	val=natural_spline(a,x,y);
	printf("Valoarea polinomului este:\n");
	disp(val);
	%solutia este:  val= 3.15


endfunction

